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We investigate the superconducting phase diagram and boundary modes for a quasi-ID system 
formed by three Fe-Chains on an s-wave superconductor, motivated by the recent Princeton ex¬ 
periment. The I • s onsite spin-orbit term, inter-chain diagonal hopping couplings, and magnetic 
disorders in the Fe-chains are shown to be crucial for the superconducting phases, which can be 
topologically trivial or nontrivial in different parameter regimes. For the topological regime a single 
Majorana and multiple Andreew bound modes are obtained in the ends of the chain, while for the 
trivial phase only low-energy Andreev bound states survive. Nontrivial symmetry reduction mech¬ 
anism induced by the I • s*term, diagonal hopping couplings, and magnetic disorder is uncovered to 
interpret the present results. Our study also implies that the zero-bias peak observed in the recent 
experiment may or may not reflect the Majorana zero modes in the end of the Fe-chains. 
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Introduction - The search for non-Abelian Majorana 
zero modes (MZMs) [IH3], which have potential applica¬ 
tions to fault-tolerant topological quantum computation 
Hi], is a focus of research in condensed matter physics. 
MZM exists in the vortex core of a two-dimensional (2D) 
(p + ip)-wave topological superconductor (SC) [8|, and 
at the end of a ID p-wave SC [9|. Theoretical propos¬ 
als showed that topological superconductivity can be ob¬ 
tained with heterostructures formed by conventional s- 
wave SC and topological insulators uni or semiconduc¬ 
tors with a Zeeman splitting [nHE]. In such devices, the 
spin-orbit (SO) interaction drives the original s-wave SC 
into an effective p-wave SC, leading to MZMs when the 
system is in topologically nontrivial regime. Motivated 
by these proposals, recent experiments using semicon¬ 
ducting nanowire/s-wave SC heterostructures observed 
the zero-bias peak (ZBP) in the differential tunneling 
conductance spectra [ 21 H 22 ], which is a suggestive sig¬ 
nature of MZMs pTH2Q] . while different theoretical in¬ 
terpretations are also available for the ZBP observed in 
the experiments [24l [25| . 

Very recently, a spatially resolved ZBP is observed by 
scanning tunneling spectroscopy (STM) in the chains of 
Fe atoms which exhibit ferromagnetic ordering and are 
placed on the surface of an 5 -wave SC (Pb) [26]. This 
study is motivated by but different from the earlier pro¬ 
posals of realizing ID topological SC by adatoms with he¬ 
lical spin configurations on an 5 -wave SC m- It has been 
interpreted in theory that such Fe-chains may exhibit 
topological superconductivity since the Zeeman splitting, 
5 -wave SC order, and Rashba SO interaction can be in¬ 
duced through the couplings between Fe atoms and Pb 
SC substrate [26l[28]. Nevertheless, the valence electrons 
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FIG. 1: (Color online) Ferromagnetic Fe atomic chains 
stacked on the Pb substrate SC along x (a) and y (b) di¬ 
rections. The later configuration may be more relevant 
for the recent experiment [26]; (c) The hopping couplings 
between nearest-neighbor, next-nearest-neighbor, and next- 
next-nearest-neighbor sites {tfjC = 1,2,3) are considered. 


of Fe atoms occupy the d-orbital states which bring about 
ten bands for a single Fe-chain. Thus for the quasi-ID 
system formed by three Fe-Chains, as considered in the 
experiment, in general there are large number of bands 
crossing the Fermi energy, which may result in compli¬ 
cated phase diagram and boundary modes, besides the 
possible MZMs suggested in the experiment. 

In this letter, we investigate the superconducting phase 
diagram and boundary modes for the three Fe-Chains 
placed on the Pb s-wave SC [Fig. (a-c)]. We find that 
the I • s term for d-orbital electrons, inter-chain diag¬ 
onal hopping couplings, and magnetic disorders in the 
Fe-chains play crucial roles in determining the symme¬ 
try classes of the superconducting phases. For the topo¬ 
logical phase we show that a single MZM and multiple 
Andreev bound states (ABSs) are obtained in the ends 
of the chain, while in the trivial regime only low-energy 
ABSs exist. The novel symmetry reduction mechanism 
is uncovered clearly to interpret these results. 
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Model- In general the configuration of the Fe atomic 
chains on the Pb surface may be complicated due to the 
strong Fe-Pb bonding. To capture the essential physics, 
we consider here the simplest situation that the three Fe 
chains sit in parallel on the surface and along ^ direction, 
with intra-chain and inter-chain hopping couplings being 
taken into account. In particular, without loss of general¬ 
ity, we consider here the hopping couplings between up to 
next-next-nearest neighbor sites [Fig. [^(c)]. Moreover, 
Two different configurations with chains stacked along 
X and y directions [Fig. (a,b)] are considered. With 
the Slater-Koster basis [29] the effective model without 
disorder can be described by 

-^triple ~ f^F EE cL(r)ca<T(r) + 

a r,cr 

+ E E E ^aa'(r>r')4<T(r)cav(r') + h.C. 

a,a' j,a r^r' 

+ E E 4<T(r)^sL'Ca<T'(r) 

a r,aa' 

TiAgo ^ ^ ^ ^ 

a,a' r,a,a' 

^ + cte^) + h.C. 

Oi r,cr,cr' 

+ E [As(r)cat(4ca^(r) + h.C.] , (1) 

a,r 

where Caa annihilation (creation) operator 

for the d-orbital electron, with a and a =t, i labeling the 
five orbital states and spin states, respectively, /ii? is the 
chemical potential, J is the Stoner-theory spin splitting 
energy with magnetization along y direction. Ago repre¬ 
sents the onsite SO coefficient, Ir is the Rashba SO coef¬ 
ficient, and Ag is the proximity induced s-wave SC order 
in the Fe chains. The coefficients (r, r') represent the 
nearest-neighbor (for j = 1), next-nearest-neighbor or di¬ 
agonal (for j = 2), and next-next-nearest-neighbor (for 
j = 3) hopping couplings in the Slater-Koster approx¬ 
imation [Fig. (b)]. The numerical magnitudes of the 
parameters are given in the Supplementary Material [39] 
according to the density function theory calculation in 
Ref. [26]. The spin splitting J ^ 2.7eV is the largest en¬ 
ergy scale in the above formula, which leads to the full 
spin polarized bands in the Fermi energy. In the presence 
of the Rashba SO interaction induced by the interface hy¬ 
bridization between Fe d-orbital and Pb p-orbital states, 
the proximity induced s-wave SC may be driven into an 
effective p-wave topological SC, which can support MZMs 
at the Fe-chain ends [IMS]. 

From the Slater-Koster basis, one can check that the 
onsite I • s and inter-chain hopping terms can 

induce the couplings between different d-orbital states. 
These couplings can qualitatively affect the symmetries 
respected by the Hamiltonian. We shall show that the 


these coupling terms and the magnetic disorder, which 
shall be considered later, play the crucial roles in deter¬ 
mining the symmetry classes and thus the superconduct¬ 
ing phases. Without such terms, the system generically 
has multiple MZMs localized in each end of Fe-chains, 
while the presence of such couplings can mix the MZMs, 
giving rise to low-energy ABSs. For the sake of a clear 
understanding of these effects, in the following we exam¬ 
ine the present Fe-chain system step by step. 

Single Fe-ehain ease.- Let us first consider the sim¬ 
plest situation with a single Fe Chain. This is equivalent 
to study the Hamiltonian 0 without inter-chain cou¬ 
plings. The single chain Hamiltonian in the momentum 
space reads Hsingie(Aso: ^z) = ( 2 Vi cos kz + 2 V 3 cosSkz — 
hF)Tz + 2tRS\iikzSx - J/2sy + AgSyTy + Ago/ • STz, where 
Vi and V 2 are matrices corresponding to hopping terms 
and respectively [39]. A key feature of the 
present system is that the symmetry class of the Hamil¬ 
tonian depends on onsite SO term. In the absence of 
I • s term, namely, if Ago = 0, we find that the Hamil¬ 
tonian respects both the time-reversal (TR) symme¬ 
try T and charge conjugation symmetry C defined via 
FT^single(Aso ~ 0? ^2;)^ = 7 ^single(Aso ~ O5 ^z\ and 

CHsingle(Aso = 0, kz)C~^ = -^*ingie(Aso = 0, -k^), with 

T = Ks,t„ C = T^, F=C^ = 1. (2) 

Here K is the complex conjugate operator. The above re¬ 
sult implies that Hgingie(Aso = ^^kz) belongs to the BDI 
symmetry class according to the ten-fold topological clas¬ 
sification [3Qll33] . which can protect integer number of 
MZMs at each end [34H38] [see Fig.[^(a)]. Note that the 
symmetries T, C do not transform orbital states, implying 
that each orbital band at the Fermi energy contributes 
one MZM. Thus the multiple MZMs correspond to the 
multiple orbital subbands crossing the Fermi energy. 

Once the I • s term is switched on, the TR symmetry 
T defined in Eq. 0 is broken. Indeed the multiple d- 
orbital states can be equivalently treated as the degree 
of freedom of a “synthetic” transverse dimension. Then 
the I • s and ID Rashba terms render the SO couplings 
in the synthetic and physical dimensions, respectively, 
giving an effective 2 D SO coupling which cannot be real 
and thus breaks the above T symmetry. Only the charge 
conjugation symmetry keeps and the symmetry class of 
the system is reduced from BDI class to D class, with 
the topology being classified by a E 2 invariant, calcu¬ 
lated by (-1)"" = sgn{Pf[H{kz = 0)rx]Ff[H{kz = F)rx]}. 
The topologically nontrivial (trivial) phase corresponds 
to u = 1(0). From Fig. (b) we can see that while in 
the most region the phase is topological, there are small 
regions which are topologically trivial. In the topological 
phase, only a single MZM is obtained in each end of the 
chain, with several ABSs coexisting. These ABSs origi¬ 
nate from the mixing between MZMs obtained in Fig. 
(a) by the I • s symmetry-breaking term. In the trivial 
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FIG. 2: (Color online) Phase diagram and boundary modes 
for a single Fe chain with (a) Ago = 0 and (b) Ago = 0.15eV 
for the onsite I • s term. Other parameters are taken as Ir — 
O.leV, As = 0.04eV, and the remaining hopping coefficients 
are given in Supplementary Material |39j . The trivial phase 
corresponds to tt-mf = 0 in the notation [umf ,'^^abss] • 



FIG. 3: (Color online) Phase diagrams and boundary modes 
for the triple Fe-chain case with (a) = 0 and (b) = 

(and disorder amplitude (^J^ax = lOmeV for the con¬ 
figuration (b) of Fig.0. Other parameters are Ago = GOmeV, 
tn = O.leV, As = 0.04eV, and the remaining hopping coeffi¬ 
cients are given in Supplementary Material [39] . 


phase, only low-energy ABSs are obtained. 

Triple Fe-chain case - Now we turn to the three Fe- 
chain model given in Eq. 0- We first consider the config¬ 
uration (a). It will be shown that in this case the diagonal 
hopping coupling term also becomes crucial in determin¬ 
ing the phases. To see this effect clearly, we parameter- 
ize the three-chain Bloch Hamiltonian Htripie(^^ao 
function of It is interesting that a new set of TR 

(T) and charge conjugation (C) symmetries are found 
when = 0, satisfying T'Htripie(^ai/ = 0,kz)T~'^ = 
^tripie(iai/ = 0,and cnuipieital = 0,kz)C~^ = 
“^tripie(*La' = 0’ ~kz)- Here C = C, while the TR sym- 
metry is given by 


f = UU^KszTz, 



z 0 0 0 

0 i 0 —i 

0 0 v/2 0 

0 10 1 

10 0 0 


—i 

0 

0 

0 

1 


( 3 ) 


Here is a local unitary matrix acting on the five d- 
orbital bases [dxy^ dxz^d^^^dyz, and one can ver¬ 
ify that = 1. Thus the Hamiltonian ~ 

belongs to a new BDI symmetry class characterized by 
T and C, and can support multiple MZMs. On the other 
hand, since the onsite I • 5 term breaks the TR symmetry 
T as defined in the single-chain model and leads to low- 
energy ABSs, in general there are both multiple MZMs 
and multiple ABSs in the present three-chain system, as 
shown numerically in Fig. (a). The ABSs are due to 
the couplings in MZMs induced by the I • s SO term. 

Similar as the result in the single-chain system, the di¬ 


agonal hopping in configuration (a) can break the new 
TR symmetry of Htripie(^aa' = 0,/c^) by verifying that 

^’^triple(^i2/5 ^ 7^ ^triple(^i2/5—fe) fol* ^ 

This implies that the symmetry class of the system is 
again reduced from the new BDI class to the D class, 
and the topology is classified by the Z 2 invariant. There¬ 
fore, the diagonal hopping terms can couple the remain¬ 
ing multiple MZMs, with only a single MZM surviving if 
the total number of MZMs is odd for the regime without 
diagonal hopping couplings. On the contrary, if the num¬ 
ber of MZMs for Htripie(^i^2' ~ O 5 fe) is even, the diagonal 
hopping term drives the system into a trivial phase with 
only low-energy ABSs existing. We plot the phase dia¬ 
gram by numerical calculation in Fig. (b). It is clear 
that both topologically nontrivial and trivial phases are 
obtained in the large ranges of chemical potential, with 
each gap closing point separating a topological phase and 
a trivial phase. Due to multiple subbands crossing the 
Fermi energy, there are multiple ABSs obtained in almost 
all the different parameter regimes. 

We emphasize that for triple-chain case the symmetry¬ 
breaking mechanism also requires the magnetization to 
be perpendicular to the Fe-chain stacking direction m- 
This implies that for the second configuration [Fig.[^(b)], 
even the diagonal hopping term cannot break the T sym¬ 
metry since the averaging magnetization is parallel to the 
stacking {y) direction. Nevertheless, in the realistic sys¬ 
tem, this symmetry can be broken when random mag¬ 
netic disorder with magnetization 5J{y) = dJ^Cx along 
X direction is present, giving the disorder Hamiltonian 
Vdis = EaEr.<T<T'4a(r)(5J"^(r)s®<,,]w(r). The further 
inclusion of the magnetic disorder Vdis in the configura- 
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tion of Fig. (b) leads to the same phase diagram shown 
in Fig. [^(b). 

Minigap and tunneling spectra - With the existence of 
multiple ABSs in the end, it is important to calculate 
the minigap, defined as the energy of the lowest ABSs 
^mini = minjF'ABs}, of the real system. A sizable mini¬ 
gap is necessary to distinguish the topologically nontriv¬ 
ial phase from the trivial phase by STM measurement. 
The eigenvalues of the boundary modes at the end can 
be calculated with surface Green’s function through iter¬ 
ation methods HOI. The calculation can be performed 
by transfer matrix method. For the present ID sys¬ 
tem with next-next-nearest-neighbor hopping couplings, 
one can separate the Fe-Chains into many principle seg¬ 
ments (PSs) along ^ axis, with each PS containing q>2 
Fe atoms along the chain direction (the total Fe atom 
number in a PS is then 3q). For the case with mag¬ 
netic disorder, one needs to take 3g ^ 1 to avoid the 
numerical error [39]. The Green’s function of the sys¬ 
tem is denoted by (cj), where n, n' are the PS in¬ 
dices and /,/' (= l,...,g) denote the q atomic layers in 
each PS. The surface Green’s function corresponds to 
n = n' = / = /' = 0 , and can be solved through 


Goo(<^) 


/ 

cj + i5^ — Hqq — Hq\T 


(4) 


Here / is a unit matrix, i^oo is the block Hamiltonian of 
the surface PS, and i^oi represents the couplings between 
the surface PS and the next PS which include the hopping 
couplings and the Rashba SO term. The transfer matrix 
T is obtained by iteration method [39] 

Nc n—1 

T{uj) = to + '^ tjrn, (5) 

n=l m=0 

where to = {uj - Hoo)~^hI^Jo = (cj - CLrn = 

(/ tjYi—\ijYi—\ ijYi—\tjYi—\) a^_^ with a^ji = tjyi or ijyi^ 
and Nc is a cut-off. The local density of states are then 
obtained by 


n, 1) =-^[G^|^(u;)], n = / = 0 , 

TT 


(6) 


with which one can determine the spectra of the bound¬ 
ary modes. With the obtained surface Green’s function 
one can further calculate the tunneling conductance by 
considering the tunneling coupling between a metallic 
lead and the end area of the chains. 

The numerical results in the experimentally relevant 
parameter regimes are shown in Fig. where (a) and 
(b) shows the minigap as functions of the the ratio r] = 
between the diagonal hopping and nearest- 
neighbor hopping strengths with fixed t ^^^,, and the am¬ 
plitude of the random magnetic disorder respec¬ 

tively. The physical regime corresponds to 77 = 1 under 
the Slater-Koster approximation m- The coefficient of 



E/Ebuik Ej Ebuik 


FIG. 4: (Color online) Lowest ABSs for the trivial phase with 
Ago = 60meV and p — — l.OeV. Minigap as a function of g (a), 
and disorder amplitude (^J^ax with 77 = 1 (b). (c) The local 
density of states (LDOS) contributed from the lowest ABSs 
for 77 = 1 , and (d) the spatially resolved tunneling spectra 
with a tunneling energy of l.OFmini- 


the onsite SO coupling /-s is fixed with Ago = 60meV [26] . 
which may determine the minigap of the system if 77 (and 
is large. Fig. (a,b) shows that the minigap 
in a topologically trivial regime increases with the di¬ 
agonal hopping coupling, and is typically a few peV for 
77 ^ 1 and < 30meV. The local density of states 

are shown in Fig.[^(c), which reflects that no MZMs but 
low-energy ABSs exist in the topologically trivial phase. 
However, from Fig. (d) we can see that when the tun¬ 
neling energy is over the minigap, the tunneling spectra 
manifests a ZBP, which brings about challenges to dis¬ 
tinguish the topological phase from trivial phase by STM 
measurement. The more precise measurement is desired 
to identify the MZMs. 

Discussion and conclusion - Several issues are worth¬ 
while to be mentioned. Firstly, the configurations we con¬ 
sider are straight and parallel Fe-chains within a single 
{x — z 01 y — z) plane. In the real system the chains may 
be staggered or not in-plane due to complexity of the Fe- 
Pb interface. A staggered and not-in-plane configuration 
can generically bring about additional hopping couplings 
which also contribute to the symmetry reduction from 
the second type of BDI class to D class. Secondly, it was 
shown that the strong coupling between the Fe atoms 
and the substrate Pb atoms can strongly renormalize the 
Fermi velocities of the Fe bands m- This effect is essen¬ 
tially because the transverse wave functions (penetrat¬ 
ing into the substrate SG) of Fe electrons have strong 
dependence on the momentum kz due to the interface 
couplings, which greatly suppresses the dependence of 
energies on kz. Nonetheless, if the interface couplings be¬ 
tween Fe and Pb atoms do not break the TR symmetries 






















5 


T and f defined in the present work, the main results pre¬ 
dicted here shall not be affected by the renormalization 
of Fermi velocities. A more detailed study of this effect 
will be performed in the next work. Finally, The long- 
range hopping couplings along the chain and transverse 
direction do not affect the TR symmetries. Therefore, 
taking into account hopping couplings between Fe atoms 
with even larger distance does not change our results. 

In conclusion, we have studied the superconducting 
phase diagram and boundary modes for a quasi-ID sys¬ 
tem formed by three Fe-Chains on an s-wave SC. We 
uncovered a nontrivial symmetry reduction mechanism 
with two different types of BDI classes (characterized by 
Z invariant) reduced to D classes (with Z 2 invariant) by 
the onsite I • 5 term for d-orbital electrons, inter-chain di¬ 
agonal hopping, and magnetic disorder couplings, which 
governs the properties of boundary modes in the topo¬ 
logically nontrivial or trivial phases. For the topological 
regime a single Majorana and multiple Andreew bound 
modes are obtained in each end of the chain, while for 
the trivial phase only low-energy Andreev bound states 
survive. Our results call for further experimental studies 
to identify the MZMs out of low-energy ABSs. 

We thank Fa Wang, Xin Liu, and Ting-Fung Jeffrey 
Poon for helpful discussions. This work is supported in 
part by the Thousand-Young-Talent Program of China. 
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Symmetry Reduction and Boundary Modes for Fe-Chains on an s-wave Superconductor — Supplementary 

Material 


A. Single Fe-chain case 


Let the spin be polarized along y direction, the whole single chain Hamiltonian in the momentum space can be 

written as Ago, 2;) ~ "^single (^so ~ H“^so 5 where 7 ^gingie(^so ~ ^i^z) ~ (2 Vi cos 2 V3 cos 

Tz ^ So 2tR sin kzls ^ To ^ Sx — J/2 Ib ^ tq 0 Sy 2 AsIb ^Ty 0 Sy^ with Vi and V 3 are matrices corresponding to 
hopping terms and respectively, given in Tablets I and 11. is the unit matrix of band degree, r and s 


or Vi(eV) 

dxy 

dxz 

42 

dyz 

dx^ —y^ 

dxy 

-0.1445 

0 

0 

0 

0 

dxz 

0 

0.5760 

0 

0 

0 

d^2 

0 

0 

-0.6702 

0 

0 

dyz 

0 

0 

0 

0.5760 

0 

1 

to 

0 

0 

0 

0 

-0.1445 


TABLE I: Nearest-neighbor hopping couplings along z direction. 


or 2 ''y 3 (eV) 

dxy 

dxz 

d^2 

dyz 

dx"^ —y‘^ 

dxy 

-0.1445 

0 

0 

0 

0 

dxz 

0 

0.5760 

0 

0 

0 

dz2 

0 

0 

-0.6702 

0 

0 

dyz 

0 

0 

0 

0.5760 

0 

dx2 - y2 

0 

0 

0 

0 

-0.1445 


TABLE II: Next-next-nearest-neighbor hopping couplings along z direction [T]. 


represent Pauli matrices of particle-hole and spin respectively. The onsite spin-orbit (SO) coupled Hamiltonian reads 


'Hso = 


( 0 

0 

Sy 

0 

0 

0 

0 

— Sy 

— Sy 

0 

0 

0 

0 

— Sy 

0 

0 

0 

0 

^/^Sx 

0 

0 

0 

0 - 

->/3S; 

Sx 

0 

-isz 

0 

0 

Sx 

0 

-iSz 

2isz 

0 

Sx 

0 

0 

2isz 

0 

Sx 

n the Nambu 

space. 



0 

0 

y/SSx 

0 

0 

0 

-V35 

0 

0 

0 


0 

0 

0 

0 

0 

0 

V^Sy 

0 

0 


-2isz 


\ 


Note that Hso is writt 

When there is no onsite SO interaction, the Hamiltonian 'Hsingie(As, 
symmetry T and charge conjugation symmetry C defined via THsingie(^so = ^ikz)T~^ = Hsingie(^so = 0? ~kz)i and 
CHsingle(Aso = 0 , kz)C-'^ = -^single(Aso = 0 , -kz), with 


T = KIb 0rz0 Sz, C = Ib 0rx0 sq , T 


0 

Sx 

0 

—2is. 

iSz 

0 

Sx 

0 

0 

iSz 

0 

Sx 

V^Sy 

, 0 

0 

0 

0 

's/^Sy 

0 

0 

0 

0 

Sy 

0 

0 

0 

0 

Sy 

Sy 

0 

0 

0 

0 

Sy 

0 

0 

= ^^kz) respects 

both the til 

^so 

0,fc.)T-i 

^single (As 

’ = c 

2 = 1. 






{SI) 


(S2) 


Here K is the complex conjugate operator. The above result implies that Hsingie(Aso = ^^kz) belongs to the BDI 
symmetry class. 

It is easy to check that T'HsoT~^ 7 ^ Hso, Cl-LsoC~^ = —Hgo, thus the TR symmetry T defined in Eq.(2) is broken 
by the onsite I • s term. Only the charge conjugation symmetry keeps. As a result, when I • s term is present, the 
symmetry class of the system is reduced from BDI class to D class. 
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B. Triple Fe-chain for the configuration (a) 


For the configuration (a) the spin is polarized along y direction, and the chains are stacked along x direction. We 

have T'Htriple(^at = 0,kz)T~'^ = = 0,-kz), and C'Htriple(^at = 0,kz)C~'^ 

where 


-^triple(C = 0,-fc.), 


f = UKUhzTz, 


i 0 0 0 —i 

0 i 0 —iO 
0 0 72 0 0 
0 10 10 
1 0 0 0 1 


(S3) 


Here U is a local unitary matrix acting on the five SK d-orbital bases [dxy^dxz^dz‘^^dyz^dx‘ 2 -y 2 ]^. The physical 
meaning of U is that it transforms between the SK bases and the bases [+2,+1,0, —1, —2]^, together with a 

local 7 r/ 2 -rotation on the d-orbital states with respect to 2 ; axis. The whole symmetry operator T includes an onsite 
spatial reflection along x direction. The diagonal hopping matrix is given by Tablet III [T] From this tablet one 


(v^)^+(eV) 

dxy 

dxz 

dz2 

dyz 

dx2 y2 

dxy 

0.2158 

0 

0 

-0.3603 

0 

dxz 

0 

-0.5388 

0.1138 

0 

0.1971 


0 

0.1138 

0.3630 

0 

-0.3689 

dyz 

-0.3603 

0 

0 

0.2158 

0 

dx2 - y2 

0 

0.1971 

-0.3689 

0 

-0.0629 


TABLE III: Diagonal hopping couplings for configuration (a) [T]. Here j i + acx + acy. 


can see that the diagonal hopping terms mix dxy^dxz states and other states, and such mixing explicitly breaks the 
aforementioned spatial reflection. One can verify that ^ 7 with 7 O 7 

the T symmetry is broken. The new BDI class is reduced to D class by the diagonal hopping couplings. 


C. Triple Fe-chain for the configuration (b) 


For the configuration (b) the spin is polarized along y direction, and the chains are also stacked along y direction. 
In this configuration the magnetization is in-plane, and the diagonal hopping matrix is given in Tablet IV [T] The 


(v^)®+(eV) 

I'J 

dxy 

dxz 

dz2 

dyz 

to 

1 

to 

dxy 

0.2158 

0.3603 

0 

0 

0 

dxz 

0.3603 

0.2158 

0 

0 

0 

dz2 

0 

0 

0.3630 

-0.1138 

0.3689 

dyz 

0 

0 

-0.1138 

-0.5388 

0.1971 

dx2 - y2 

0 

0 

0.3689 

0.1971 

-0.0629 


TABLE IV: Diagonal hopping couplings for configuration (b) [T]. Here j i + acx + aey. 

hopping terms are block diagonal and do not mix dxy^dxz states and other other states. Thus the diagonal hopping 
couplings cannot break the onsite spatial reflection along x direction, and accordingly do not break the T symmetry. 
On the other hand, the nearest- and next-next-nearest-neighbor couplings along y direction take the forms in Tablets 
V and VI. These terms also do not break the T symmetry. 

In the realistic system, a random magnetic disorder can exist with nonzero magnetization 5J{y) = dJ^ix along x 
direction, which gives the disorder Hamiltonian 

a r,crcr' 


(S4) 


























(eV) 

1.7 

dxy 

dxz 

d^2 

dyz 

1 

to 

dxy 

0.5760 

0 

0 

0 

0 

dxz 

0 

-0.1445 

0 

0 

0 

d^2 

0 

0 

-0.2759 

0 

-0.2276 

dyz 

0 

0 

0 

0.5760 

0 

dx‘2 - y2 

0 

0 

-0.2276 

0 

-0.5388 


TABLE V: Nearest-neighbor hopping coupling in y direction. Here j = i acy. 


(eV) 

dxy 

dxz 

dz2 

dyz 

dx^ — 

dxy 

0.5760 

0 

0 

0 

0 

dxz 

0 

-0.1445 

0 

0 

0 

dz2 

0 

0 

-0.2759 

0 

-0.2276 

dyz 

0 

0 

0 

0.5760 

0 

dx2 - y2 

0 

0 

-0.2276 

0 

-0.5388 


TABLE VI: Next-next-nearest-neighbor hopping coupling in y direction [T]. Here j — i -\- 2aey. 


This term transforms according to TVdis^~^ = — Vdis- Therefore, the further inclusion of the magnetic disorder Vdis 
breaks the T symmetry in the configuration (b). Note that the disorder with in-plane magnetization (along y and ^ 
directions) does not affect the symmetry. 


[1] W. A. Harrison, Electronic Structure and the properties of solid: The Physics of the Chemical Bond, Dover Publications 

























